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The background field method has been used successfully to determine hadron electromagnetic 
polarizabilities. Recently questions have been raised regarding the proper way to deal with the 
electric field on the lattice. In this paper, we show how the presence of a background electric field 
affects the quark hopping matrix. We use this formulation to carry out simulations on quenched 
configurations, and we present the results for neutron electric polarizability for pion masses as 
low as 500 MeV. We find that the polarizability is roughly constant for the quark masses consid- 
ered, a ~ 1.5 X 10^'* fm^. While the polarizability is positive, it is significantly smaller than the 
experimental value due to the fact that the quark masses used are too large. 
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1. Introduction 

Background field method has shown promising results in computing hadron electromagnetic 
moments and polarizabilities [1, 2, 3, 4, 5, 6, 7]. The basic idea [1] is to introduce a static electro- 
magnetic field on the lattice and then measure the change in the hadron masses due to the presence 
of the field. To the lowest orders, the shift in the hadron masses is given by 

^E = -p-E-'jx-B-]^{aE'^ + fiB^) + ..., (1.1) 

where p and Ji are the electric and magnetic dipole moments and a and j3 are the electric and 
magnetic polarizabilities. To introduce the field, the charged particles get coupled minimally with 
the electromagnetic field; the covariant derivative becomes 

D^ = d^- igGpL -iqAfx, (1.2) 

where is the chromoelectric field and is the static electromagnetic background. On the 
lattice, the presence of the electromagnetic background modifies the fermion hopping matrix; from 
a practical point of view, the change can be viewed as an additional phase factor multiplying the 
links: 

U^^e-^"''^i^U^. (1.3) 

A complication arises due to the Euclidean nature of lattice formulation; formally, one has the 
following rules for converting the Minkowski formulation to Euclidean: 

■^1,2,3^^1,2,3 Ai^2,3 ^^1,2,3 

Xq ^ X^ = iXQ Aq — > A4 = —iAq. (1-4) 

To introduce a constant electric field in the x-direction we can choose the following potential: 

Am = {0,+Et, 0,0) ^ AE = {-iEx4,0,0,0) or (1.5) 
Am = (-£x, 0,0,0) ^ Ae = {0,0,0, +iExi) , (1.6) 

where Am stands for Minkowski and Ae for Euclidean formulation; note that in the Minkowski 
formulation the components are {Ao,A) and in the Euclidean one they are (A,A4). Similarly, for a 
magnetic field in the x-direction we can use the following potential choices: 

Am = (0,0,+Sz,0) ^ A£ = (0,+Sz,0,0) or (1.7) 
Am = (0,0,0,-63;) ^ As = (0,0,-63;, 0). (1.8) 

Thus, on the lattice we can use the following phase factors to produce a constant electric or mag- 
netic field: 

E^: Ux ^ e-''"'^'''Ui or ^ e^'^'^'^'VA (1.9) 
5,: U2^e-'''"'^'''U2 or U3 ^ e'''"'^''''U3. (1.10) 

The surprising result in Eq. (1.9) is that the phase factors are real, in contrast to the magnetic case 
where the phase factors are U{1) phases. This can be traced back to the extra / factor that appear- 
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when performing the rotation to the "imaginary" time: for the first choice of potential it arises from 
the fact that the time becomes imaginary, whereas for the second choice it comes from changing Aq 
to —iA^. This result is at variance with the standard treatment [1, 2, 3, 4, 5] where both the electric 
and magnetic field are introduced using a i7(l) phase. The fact that the electric field defined using 
SiU{\) phase corresponds to an imaginary electric field in Minkowski space was first pointed out 
by E. Shintani et all [8]. 

The purpose of this paper is to show that the result presented above, derived using/onna/ rules, 
is correct. In section 2 we first argue for its correctness using the Wilson loop as a simple example 
and then carry out an explicit calculation for charged bosons placed in a constant electric field. 

In section 3 we address the issue of the exponential versus linear phase factor. It is argued [1] 
that the phase factors in Eq. (1.9) should be replaced with a linearized version. We show that it is 
better to keep them in the exponential form. 

Finally, in section 4 we present the results of our simulations for the electric polarizabilities 
of the neutron. We show how to compute its polarizability using either the real phase factor or the 
U{\) phase. We also present the result for the linear case and compare it with the exponential case. 

2. Euclidean formulation 

To show that the use of a real phase factor in the electric field case is justified imagine that we 
place a heavy quark-anti-quark pair in a electric field (see Fig. 1). For large T the Wilson loop 
average decays exponentially, i.e. (W) ~ e^^'^^^. In the presence of an electric field the energy of 
the state is altered by the dipole energy: V{R) — > V{R) — qER, where q is the charge of the quark. 
This defines our expectations. 
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Figure 1: Heavy quark-anti -quark pair in the presence of an electric field. 

When the electric field field is introduced on the lattice using the first choice in Eq. (1.9) the 
links in the x direction get an addition phase factor; the Wilson loop is then changed: the bottom 
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links in Fig. 1 contribute a factor (^e^viEt-^R/a _ ^-qERt jj^g contribute = 

^qER{t+T) _ Since this factor is the same for every configuration, the average value of the Wilson loop 
becomes 



(W) ^ {W)e 



-qERt^qER{t+T) ^ ^-V{R)T+qERT 



(2.1) 



We see then that the potential is modified according to our expectations; the change in the energy 
due to the presence of the electric field Am = —qER is exactly the dipole energy. 

To further confirm that this is the correct way to introduce the electric field we can compute 
the partition function in the Euclidean formalism for a charged bosonic field. The Lagrangian of 
the system is 



d^x [d^(j)*d^^-m^^*^] . 



(2.2) 



We introduce the background field using the minimal substitution d^^d^ + iqA^ . The Lagrangian 
in the presence of the field becomes: 



L = Jd^x [{d^-iqA^)(p*{d^+iqA^)^-m^(j)*(l)], 



(2.3) 



where we see that the covariant derivative acting on <p* is — iqA^ due to the fact that the conjugate 
field has opposite charge (it is also needed to preserve the reality of the Lagrangian). 

To get to the Euclidean formulation we need the Hamiltonian of the system which we get using 
the Legendre transform 

H = I d\[n{x)^{x) + 7i*{x)^*{x)]-L, (2.4) 
— ^*{x) — iqAo{x)^*{x) K The Hamiltonian of this system is 



where 7i(x) 



H 



d^x 



71*71 + /^Ao( 71*0* - 710) + (V-/^A)0*(V + /^A)0 +m^0*0 



(2.5) 



To quantify this system we discretize it, introduce field operators defined on lattice points satisfying 
the canonical commutation relations and then normal order it. Since we are interested in the electric 
field we will set A = 0. The resulting hamiltonian is: 



^ = L [K^n + iqAoift:^: - A,4n) + (V0„)* + m^;Jn] , 

n 

where [0„, ?!„,] = /5„,„, and V is some discretization of the gradient. 

To determine the Euclidean action we need to compute the matrix element 



(2.6) 



0r 



(2.7) 



where |0) is a vector in the "field" representation (i.e. 0„ |0) = 0„ |0) and 0* |0) = 0* |0)) and 
is the Euclidean Lagrangian density. The calculation is straightforward; the most difficult steps are 

1 -J-irf,'-, 







TT e 'f 



and 



,-a,iqAo{n*^;-1t„^„) i 



(2.8) 
(2.9) 



e "'■■2**"\"'i'^n ""I-"' 10) = |0) , 

'Note that while this is the standard representation for the charged Klein-Gordon field, the use of complex variables 
is formal; the discussion is correctly carried out in terms of two real fields 2 which are related to our field by = 



((^l+/02)/v/2. 
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where 0„, = 0„, for m^n, 0„ = e"'^^'^^n and = e "''i^o^* 2 Using the relations above we get 



The important thing to notice is the electric field contribution appearing in the hopping term in the 
time direction. There are no links in this formulation, but we see that the field propagating forward 
in time gets a factor due to the electric field e"'''"^" ; this is exactly the same factor as the one derived 
using the formal rules for the second choice of potential in Eq. (1.9). 

We have shown that in the case of charged bosons the phase factor due to an external electric 
field is real. We see no reason why this should be different for fermions so we conclude that the 
rules derived formally in the first section are correct. While we didn't show it, the hopping matrix 
for a charge boson in the presence of a magnetic field has a. U{1) factor; this is also in perfect 
agreement with the factors derived using the formal rules. 

3. Exponential vs. linear phase factor 

As mentioned in the Introduction, most of the studies carried out this far [1, 2, 3, 4] use a 
linearized form of the phase factor. It was suggested in [1] that the exponential factor e^'^"^>^ 
should be replaced by its linearized version 1 — iqaA^ ; the authors argued that the the linearized 
factor makes the Dirac matrix resemble the continuum like covariant derivative where the coupling 
with the electromagnetic field is linear. 

The difference between the linear and exponential form is of the order so it seems like 
that shouldn't make much of a difference in the continuum limit. However, the change to linear 
form also changes the action at the order E'^; the polarizability is derived from the response of the 
observables at the order E'^ thus its value is affected. In the next section we will show that, at least 
when a is around 0.1 fm, this change can be quite significant. 

It may be that in the continuum limit the polarizability is the same irrespective of how we in- 
troduce the electric field. However, this is not guaranteed; for example, when computing the energy 
density of a free fermionic gas at non-zero chemical potential the lattice result is wrong when the 
chemical potential is introduced using an expression similar to the one in the continuum [9]. While 
this example might look off-topic, it is in fact very pertinent to our discussion: the chemical poten- 
tial plays a very similar role to Aq. Hasenfratz and Karsch [9] point out this similarity and argue 
that it is the gauge symmetry in Aq that protects the continuum Uke formulation from divergencies. 
They conclude that in order to avoid divergencies on the lattice, you have to couple the chemical 
potential in a gauge invariant way: the phase factors are then changed from 1 ita/x to e^"^. 

While it is not entirely clear that their argument can be carried over to our situation, it is worth 
pointing out that when using the linearized version the gauge symmetry in the external field is 
lost. To be more precise, take the two choices for the electromagnetic potential in Eqs. (1.5) and 
(1.6); they both represent an electric field in the x-direction and we expect that the results that they 
produce are the same. If we are to use the con^esponding exponential factors in Eq. (1.9) then the 

^Note that ^* is not tlie complex conjugate of this confusion is due to our improper use of complex variables. 
We treat ^ and (^* as independent variables. 



at 



(2.10) 
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results come out to be the same since the gauge symmetry is preserved. If we choose to linearize 
these factors then the results of our simulations would be different. The Wilson loop in the presence 
of the electric field introduced using the linear factors is: 



depending on whether we introduced it using the space links U\ or the time links Ui\. Notice that 
when we use the electric field using time links the Wilson loop decays exponentially in time and the 
only difference is an a dependent contribution to the dipole energy; more worrisome is the fact that 
when using the space links the Wilson loop doesn't decay exponentially with time but it acquires 
a piece that goes like . This would make it impossible to determine the energy of the state by 
simply fitting for an exponential behavior. 

We conclude that it is better to use the exponential form rather than the linear form. Even if 
the continuum limit is the same for both formulations (which is not guaranteed), the exponential 
form has the benefit of being gauge invariant in the external field. 

4. Numerical results 

To check these ideas we ran a set of quenched simulations on 24^ lattices. We used Wilson 
action at j3 = 6.0 which corresponds to a lattice spacing of a = 0.093 fm. We used clover fermions 
and our lowest pion mass was 500 MeV. For the electric field we used T] = a^qE = 0.00576, where 
q is the charge of the down quark. To avoid dealing with field discontinuity we used Dirichlet 
boundary conditions in both time and x direction. 

We computed the correlation function for the neutron GE{t) in the presence of the electric field 
E. To extract the polarizability we computed the ratio GE{t)/Go{t) which is expected to behave 
asymptotically as e^^"". Unfortunately, our lattice was not long enough in the time direction and 
we used a two exponential form to fit Ge and Gq. 

After extracting the mass shift we computed the polarizability. In Fig. 2 we plot the polariz- 
ability as a function of the pion mass. Note that the polarizability is roughly constant in the range 
of quark masses we studied. Also plotted is a result that includes dynamical fermion effects and the 
effect of fermions on the background [5] ^. We see that our result is in good agreement with this 
value and that the dynamical effects do not play an important role when the pion mass is around 
650 MeV; this is going to change when we approach the chiral limit. 

We also run a set of simulations using the same parameters but using a.U{l) factor to introduce 
the electric field - this corresponds to introducing a imaginary electric field iE. The expectation was 
that the mass shift is going to change sign since it depends on the square of the electric field. We 
verified that this is indeed true and the polarization derived using these mass shift agreed perfectly 
with the one derived using the real phase factor. 

The final set of simulations used a linear expression for the electric field contribution. We 
argued in the previous section against using it, but we wanted to check whether the difference is 

^The value reported in [5] is the negative of tiie value plotted here - this is due to the confusion regarding the proper 
way of introducing the electric field in Euclidean formulation. 



Ui : (W) 
U4 : (W) 



V{R)T+qERT^^_aq^E^R[{t+Tf-+t-] 
V{R)T+qERT ^^aq^E~T[{x+R)^+x^] 



(3.1) 
(3.2) 
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Neutron polarizability vs lUj, 
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Figure 2: Neutron polarizability as a function of the pion mass. The result in red includes the dynamical 
effects of the fermions and the influence of the electric field on the vacuum [5]. 
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Figure 3: Neutron polarizability as a function of the pion mass when introducing the electric field via a 
linear factor. 



significant. In Fig. 3 we plot the polarizability computed with the linear field. We first note that 
these values agree well with the results reported in [2] except for the sign; the sign difference arises 
from the confusion regarding the electric field in Euclidean time. We see that these results no longer 
look like they are converging towards the experimental value - even the sign differs. We also see 
that the difference between the linear and exponential case is quite significant; in the unlikely case 
that they converge to the same value it would take extremely fine lattices given the discrepancy for 
a = 0.093fm. 

Our numerical simulations support the conclusion of our previous section: to introduce the 
electric field we need to use an exponential phase factor; the linear case seems to differ substantially 
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from the exponential case. It is also important to note that we can compute the polarizability using 
a ?7(1) factor as long as we remember that this corresponds to an imaginary electric field. 

5. Conclusion and outlook 

Our goal was to understand how to relate the electric field in the Minkowski space - the 
physical electric field - to the phase factors that appear in the hopping matrix. We have shown 
that the presence of an external electric field gives rise to an exponential phase factor e^^^^' . To 
compute the polarizability we can use either the real phase factors or U{\) ones as long as we 
correct for the sign: 



The presence of the magnetic field gives rise to a ?7(1) factor and to compute the magnetic polar- 
izability we use 



Our numerical results agree well with previous simulations if we account for the sign change. 

The nucleon polarizability for m,i = SOOMeV is positive but significantly smaller than the exper- 
imental value. This is actually in agreement with chiral perturbation theory that predicts a. l/niji 
raise in polarizability as one approaches the chiral limit. We don't see the trend yet and most likely 
we need to carry out simulations at smaller quark masses. The other major challenge is to include 
the contribution of the electric field on the vacuum which we expect to have an even bigger impact 
on the value of the polarizability than the inclusion of dynamical fermions. Our current plan is to 
include this effect using a re- weighting of the determinant. 
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